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Abstract 

We study the precession of an off-axis straight vortex in a rotating non- 
axisymmetric harmonic trap in the Thomas-Fermi (TF) regime. A time- 
dependent variational Lagrangian analysis yields the dynamical equations of 
the vortex and the precessional angular velocity in two-dimensional (2D) and 
three-dimensional (3D) condensates. 
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I. INTRODUCTION 



Since the original experimental observation of Bose-Einstein condensation in dilute 
atomic gases the creation and detection of vortices has attracted great interest (for a 

recent review article, see, for example Ref. 0]). In particular, one experiment has measured 
the precession frequency uj p of a vortex line in a nonrotating stationary spherical condensate 
|| and compared u p to various theoretical estimates @-ll|. Although these theoretical ap- 
proaches largely predict the same results, perhaps the most transparent and physical picture 
||7]|| relies on a variational Lagrangian method |12| that focuses directly on the position of 
the vortex core. For an axisymmetric condensate in rotational equilibrium at a small angular 
velocity fi <Cwi, the original precession frequency u p is altered to uj p (Q) = u p — fl by the 
external rotation. The present work describes the nontrivial generalization of the variational 
Lagrangian method to a large rotating nonaxisymmetric condensate in the Thomas-Fermi 
(TF) limit. For extreme asymmetry, we may note that the precession frequency uj p (Q) be- 
comes nearly independent of Q (for Q <C ou±) because the irrotational flow induced by the 
rotating asymmetry screens the effect of the external rotation. 



II. PRECESSIONAL DYNAMICS OF A STRAIGHT VORTEX 

A Bose condensate is described by a macroscopic order parameter (the condensate wave 
function) with particle density n = \^\ 2 - Consider a condensate in a nonaxisymmetric 
harmonic trap potential Vt r (r) = ^M{oj 2 x 2 + uj 2 y 2 + uj 2 z 2 ). The condensate experiences 
both the trap potential V tT and the self-consistent Hartree interaction Vh = g\^\ 2 = gn 
arising from the interaction with all the other particles, where g = Airh 2 a/M is the effective 
interparticle interaction strength and a > is the s-wave scattering length []I3[ , |I4| . At zero 
temperature, the order parameter obeys the time-dependent Gross-Pitaevskii (GP) equation 

When the trap rotates with angular velocity Q about the z axis, the corresponding GP 
equation in a co-rotating frame acquires an additional term —QL z ty = ihQr x V\& on the 
right-hand side, arising from the transformation to the rotating frame JT3 . 



Instead of working directly with the time-dependent GP equation (|l|), it is often prefer- 
able to introduce a Lagrangian formalism based on the Lagrangian functional 

C[^] = T[^]-J=\% (2) 

where 

is the time-dependent part of the Lagrangian that is analogous to the kinetic energy in 
classical particle mechanics. Similarly, the free energy J 7 ^} here plays the role of the 
potential energy. In a nonrotating system, JF is simply the energy functional £ , which 
includes the gradient (bending) energy, the trap energy, and the self-consistent Hartree 
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energy. When the system rotates, however, the free energy has an extra term JF = 8 



flL z |15[, and the appropriate free-energy functional of the system is 

F[9] = Jdv(^\V^\ 2 + V t M 2 + ^m 4 + ^^^y (4) 

Given the Lagrangian, the associated action is the time integral S[$] = J t * 2 dtC[^}, and it is 
straightforward to verify that the action is stationary with respect to small variations of ^ 
and \&* when ^ obeys the time-dependent GP equation ([!]). Thus the time- dependent GP 
equation is the Euler-Lagrange equation for this problem. 

This Lagrangian formalism is exact, but it also provides the basis for a powerful ap- 
proximate variational method. If the wave function \1/ depends on various parameters, the 
integral in Eq. (|3]) will involve the first time derivatives of these parameters. Consequently, 
the resulting £ will serve as the effective Lagrangian that provides the corresponding dynam- 
ical equations for these parameters. This method has been used to study the quadrupole 
oscillations of a vortex- free condensate [|1^]. Here, we use it to study the dynamics of a 



vortex in a rotating nonaxisymmetric trap. 

Specifically, we assume that a straight singly quantized vortex is displaced from the 
center of the trap with transverse coordinates Xq and j/ j these variables will serve as time- 
dependent variational parameters that obey Lagrange's equations [0)0- in the TF limit, 



the vortex induces negligible change in the condensate density, and, in the presence of the 
rotation, |\1/| 2 is given by P^JT7| 



g\V(r)\ 2 = ji- V tT (v) =p,-\M (u 2 x x 2 + u 2 y y 2 + u 2 z z 2 ) (5) 

where uj x and u y are effective oscillator frequencies that determine the shape of the conden- 
sate in the rotating asymmetric trap. Specifically, they are given in terms of a parameter 
a: 

u 2 x = u 2 x + a 2 - 2aQ, Co 2 y = u? y + a 2 + 2aQ, (6) 
where a satisfies a cubic equation 

2a 3 + a(u£ + «J - Atf) + - u 2 y ) = 0. (7) 



The parameter a has a simple physical interpretation [[L6] , for it determines the irrotational 
flow vo induced by the rotating nonaxisymmetric trap through the velocity potential $o — 
axy, with v = V$ - As noted below, a = — (u) 2 — u 2 ) / \ Cj 2 + uj 2 ^j. 

Equation ([5]) shows that the condensate density has the familiar TF parabolic form 

|*| 2 = n (l-£-£-£V (8) 
°{ Rl R 2 R 2 ) 

if the right-hand side is positive and zero otherwise. Here no = jx/g = Mjl/Anah 2 is the 



central density, Ri = J2jl/Mojf (for % = x,y), R z = yz/z/Mu; 2 are the TF radii, and 
p, = ^hujQ (lbNa/do^J ^ with lJq = uj x uj y oj z and d\ = Ti/Mujq. For a small angular velocity 
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< ui where uj\ = + ou 2 ), the renormalization of the oscillator frequencies and TF 
radii is of order Q 2 and thus small. 

In our variational approach, we assume that the dominant vortex-induced contribution 
to the Lagrangian arises from the condensate's superffuid motion and use the following trial 
wave function 

# = |#| e i(5o+5l) , (9) 

where |^| follows from the TF density profile in Eq. (^). The first term of the phase is 
taken as So = (Ma/h) xy, which represents the irrotational motion induced by the rotating 
nonaxisymmetric trap (note that a vanishes if u x = uj y ) . The second term is taken to be 

S\ = arctan ( — — — ) , (10) 

\X — XqJ 

which characterizes the circulating flow around the vortex line. In contrast to the analysis 
of Refs. 0,0, we do not include an image vortex because the form of the TF condensate 



density ensures that the particle current vanishes automatically at the TF surface. 

Variation of the functional (|j) yields an Euler-Lagrange equation for So, whose solution 
gives 0530029 

ft Wx + ^yJ ft \R 2 x +R 2 y J 

Note that the time-dependent parameters xo and yo appear only in the phase Si, so that 

d& = i VdtSi 

= -i^ro-VSi 

(x - x ) 2 + {y- y ) 2 



III. TWO-DIMENSIONAL CONDENSATE 

For simplicity, we first consider a two-dimensional condensate that is unbounded in the z 
direction (taking u z — > 0), in which case the various terms in the Lagrangian are interpreted 
per unit length. To evaluate the time-dependent part (H), it is convenient to introduce 
dimensionless coordinates scaled with the renormalized TF radii and shift the origin of 
coordinates to the position of the vortex. In this way, we find 



T[tt] = hn R 2 x R 2 y I d 2 r -Z 



hrioR 2 R 2 



-ioV + VoX r 



R 2 x x 2 + R 2 y y 2 
-Xo sin + 2/o cos « 



1 — (r cos (ft + x ) 2 — (r sin (f> + y ) 



R 2 cos 2 



+ R 2 y sin 2 



2 , 
—B 3 
3 



AB 



(13) 



where 
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A = I — xl — i/q, 

B = x cos 4> + yo sin i 



(14) 



with x and j/ also dimensionless. As anticipated in the general discussion, T indeed depends 
on the first time derivatives of the parameters x$ and yo- 
Lagrange's equations now take the familiar form 



V x (J -T) = 0, V y (T-F) = 0, 



where 



_ d d 
V x = — —— 

dt oxq 



_d_ 

dXr 



(15) 



(16) 



and similarly for T> y . If this operator is applied to functions of the form appearing in Eq. (|H| 
it is straightforward to verify that 



V x [(-± sin + cos <j>)f(x , y , (f))} = 
In our case, 

(cos d XQ + sin c^) A — —2B, 



-y (cos <j)d X0 + sin(pd yo ) f(x ,y , 



(cos + sin <pd yo ) B 



1. 



(17) 



(18) 



Applying this operator to the time-dependent term (|13j) simplifies the integral considerably 
and yields 



V x T[y]=y Q hnoRlRl 



A 



Rl cos 2 6 + R 2 , sin 



2/o 2irhn R x R y (l-x 2 Q -yl) 



The first of Lagrange's equations (15) then becomes 

0A.F[tf] 



(19) 



(20) 



where AjF denotes the extra free energy associated with the presence of a vortex. This 
equation, along with a similar one for the yo dependence can be rewritten as 



2nTi no R x R y z x r = Vq A^ r (r ), 



(21) 



where r = [xoiVo] is the dimensionless coordinate vector of the displaced vortex. As is 
familiar in the context of two-dimensional vortices in superfluid helium |2]J , this dynamical 
equation can be interpreted as a balance between the Magnus force and the gradient of the 
free energy with respect to the coordinates of the vortex PJTyrO]]. 

We can calculate the free-energy part with logarithmic accuracy |18| by substituting 
([]) -([□]) into the energy functional (||): 



In ^ 



MQ R 2 x R 2 y 



(22) 
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where R\ = 2R 2 x R 2 y /(R 2 x + R 2 y ) and £ 2 = h 2 /2Mfi. With these results, we can derive the 
equations of motion and obtain the angular frequency u} p (Q) of the vortex precession. An 
easy calculation yields 

x = -u p [Q)y , y = +uj p {Q)x , (23) 

where 

uJQ) = u° v - 2RxR l Q (24) 
pK 1 p R 2 x + R 2 y 



and 



^T^M^It^ (25) 



MR X Ry 




is the precession frequency in a nonrotating condensate of stretched dimensions R x and R 



For a not-too-large angular velocity Vt <C lu±, note that u p ~ u; p (0). An earlier more intricate 
derivation of this result analyzed the dynamical motion of each element of the vortex core 
||. If ujpiVt) is positive, the vortex precesses in the positive (counterclockwise) sense, namely 
with the same sense as the circulating fluid around the vortex. Such behavior is indeed seen 
in the recent JILA experiments [[J for nonrotating condensates. 

For an axisymmetric condensate with R x = R y = R±, the precession frequency a> p (0) in 
the absence of rotation is shifted to u> p (Q) = oj p — Q by the applied rotation (as expected 
from the transformation to the rotating frame P2"|). For an asymmetric trap, however, the 
induced irrotational flow in Eq. (|TTD acts to screen the effect of the external rotation. In 
the extreme limit R y R x , for example, we have Up(Q) ~ uo p — (2R y /R x )Q, so that the 
precession then becomes nearly independent of Q for small Q. 

This behavior can be understood by examining the background fluid flow velocity in the 
rotating frame v = (h/M) V'Sq — Qz x r' in the absence of the vortex [p0| ; for clarity, we now 



use primes to denote the original dimensional coordinates with, for example, x' = xR x . In 
particular, v x (x',y') = 2Vty' R 2 X /{R 2 X + R^). This induced velocity appears in the dimensional 
equation of motion as 

_ _ u; o^o_ _|_ v x( x 'q, y'o) _ _^o^o_ _|_ l^lRx ^1 ^2g^ 
R x p Ry Rx '' Ry R 2 + Ry 

which readily reproduces the first of Eqs. (E3), and similarly for the second equation. 



IV. THREE-DIMENSIONAL CONDENSATE 

It is not difficult to generalize these results to a three-dimensional TF condensate. For 
example, the explicitly time-dependent part of the Lagrangian now contains an additional 
integral over the axial coordinate z: 
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T[tf ] = Un R 2 x R 2 y R z / d d r -~ 



3 -x y + 2/ x r 



i? 2 X 2 + # 2 ?/ 2 



2hn R 2 x R 2 y R : 
finoR 2 R 2 R z 



d 2 r 



-x y + y x 
dz — 



1 — (r cos + xo) 2 — (r sin + y ) 2 — -s 2 

1 — (r cos + x ) 2 — (r sin + y ) 2 — z<2 



Rlx 2 + Ry 



— xq sin + ?/o cos ( 
R 2 X cos 2 + -R 2 sin 2 



-5fiA 3/2 - 3B 3 A 1/2 + 3 (A + B' 



71 



— h arctan 
2 



where = 1 — (r cos + xo) 2 — (r sin + y ) 2 - Although this angular integral is difficult to 
evaluate directly, only the simpler quantity X^T^] is needed, and a straightforward analysis 
yields 



V X T[V] 



8nhn ~ ~ ~ 
Vo — ; ' ' ' 



R x R y R z (l 



„2n3/2 



(28) 



In addition, the free energy is given by || 

rl? 12 



, 87T^ 2 n R z . 
Aj^(r ,iZ) = (1 



" T 



4Mfi R\R 2 y 

5 h Ri + m 



(29) 



The dynamical equation again reduces to ([H]), and the only difference is that the angular 
frequency has a modified numerical coefficient 



uJ£l) = ~ — — In ( —r 

A ' 2 MR x R y \ £ 



2R x Ry ^ 



(30) 



this expression agrees with an earlier result based on the method of matched asymptotic 
expansions ||. 

Note that our analysis for a three-dimensional condensate considers only a straight vortex 
line. In the present Thomas-Fermi limit of a large condensate, the axis of the vortex line 
must lie perpendicular to the condensate's surface [§,[|. Hence the approximation of a 
straight vortex is directly applicable to the central region of a disk-shaped condensate. For 
cigar-shaped condensates, in contrast, the bending of the vortex generally plays an essential 
role. Indeed, for sufficiently elongated condensates, a straight vortex on the central axis can 
become unstable with respect to bending deformations P,[TT|. 



(27) 



V. DISCUSSION AND CONCLUSIONS 

The Lagrangian functional @ contains two parts, T[^f] with explicit time dependence 
and J 7 ^?] involving the time-independent free energy in the rotating frame. If only the 
latter quantity is considered (which is the Hamiltonian in the rotating frame), the stability 
can be inferred directly by considering how T changes for small lateral displacements of 
the vortex from the central position. Such methods have been used to analyze the onset of 
metastability in a rotating nonaxisymmetric TF condensate 
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In contrast, the full Lagrangian allows a more complete dynamical description. In the 
present case, the time dependence of the vortex position r yields equations of motion (pT[) . 
which readily give the precession frequency u> p (Q) in a rotating condensate. As mentioned 
in Sec. [J, the present description does not include an image vortex. Although such an 
image vortex is needed to satisfy the appropriate boundary condition for a condensate with 
uniform density in a rigid container |19| , the present TF condensate density vanishes at the 
TF boundary so that the superffuid current automatically vanishes there. In addition, the 
image vortex has negligible effect on the energy when calculated with logarithmic accuracy. 
As emphasized in Refs. 0,0], the vortex is unstable (metastable) whenever the gradient of 
AjF(r , fi) with respect to r in Eq. ( pT|) acts to move the vortex away from (back toward) 
the center of the trap. 

The present Lagrangian approach gives a direct physical derivation of the precession 
frequency. Although the results describe only a single straight vortex displaced from the z 
axis, the analysis is far simpler than the method of matched asymptotic expansions applied 
to the full Gross- Pitaevskii equation |8|||. For a one-component condensate, it fully justifies 
the intuitive dynamical equation fl2l|) suggested earlier by McGee and Holland |L0| . It would 
be valuable to generalize the present approach to the more interesting and challenging case 
of two components, where additional restoring (buoyancy) forces have been proposed [|TU|. 
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